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ABSTRACT 

Classical novae eject significant amounts of nuclear processed material into 
the interstellar medium. Among the isotopes synthesized during such explosions, 
two radioactive nuclei deserve a particular attention: 22 Na and 26 A1. In this 
paper, we investigate the nuclear paths leading to 22 Na and 26 A1 production 
during nova outbursts by means of an implicit, hydrodynamic code that follows 
the course of the thermonuclear runaway from the onset of accretion up to the 
ejection stage. New evolutionary sequences of ONe novae have been computed, 
using updated nuclear reaction rates relevant to 22 Na and 26 A1 production. 
Special attention is focused on the role played by nuclear uncertainties within 
the NeNa and MgAl cycles in the synthesis of such radioactive species. From 
the series of hydrodynamic models, which assume upper, recommended or lower 
estimates of the reaction rates, we derive limits on the production of both 22 Na 
and 26 A1. We outline a list of nuclear reactions which deserve new experimental 
investigations in order to reduce the wide dispersion introduced by nuclear 
uncertainties in the 22 Na and 26 A1 yields. 

Subject headings: novae, cataclysmic variables - - nuclear reactions, 
nucleosynthesis, abundances 



-3- 

1. Introduction 

The thermonuclear runaway model has been successful in reproducing the gross 
features of nova outbursts. According to this widely accepted scenario, classical novae 
are produced by thermonuclear runaways (hereafter, TNRs) that take place in the white 
dwarf component of a close binary system. The large, main sequence companion overfills 
its Roche lobe, providing matter outflows through the inner Lagrangian point that lead to 
the formation of an accretion disk around the white dwarf. A fraction of this H-rich matter 
lost by the companion ultimately ends up on top of the white dwarf, where it is gradually 
compressed as accretion goes on. The piling up of matter heats the envelope up to the point 
when ignition conditions to drive a TNR are reached. 

An extended set of hydrodynamic computations of classical nova outbursts has been 
performed during the last 25 years (see Starrfield et al. 1972, for the first hydrodynamic 
study of the TNR model, and the recent papers by Kovetz & Prialnik 1997, Starrfield et al. 
1998, Jose & Hernanz 1998, and references therein), for a wide range of white dwarf masses 
and initial luminosities, mass accretion rates and initial chemical compositions. From the 
nucleosynthesis viewpoint, these computations have been able to identify two types of nova 
outbursts, those occurring in CO or in ONe white dwarfs. The latter ones have provided 
a framework for the origin of the high concentrations of Ne and more massive isotopes 
found in the spectra of some well-observed novae, such as V693 CrA 1981, V1370 Aql 1982, 
QU Vul 1984 No. 2, V838 Her 1991 or V1974 Cyg 1992 (Livio & Truran 1994; Gehrz 
et al. 1998). Among the isotopes synthesized in these so-called Ne (or ONe) novae, two 
radioactive species have raised a particular astrophysical interest: 22 Na and 26 A1. 

The potential role of 22 Na for diagnosis of nova outbursts was first suggested by 
Clayton & Hoyle (1974). It decays to a short-lived excited state of 22 Ne (with a lifetime of 
r = 3.75 yr), which de-excites to its ground state by emitting a 7-ray photon of 1.275 MeV. 
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Through this mechanism, nearby ONe novae within a few kiloparsecs from the Sun may 
provide detectable 7-ray fluxes. Several experimental verifications of this 7-ray emission 
at 1.275 MeV from nearby novae have been attempted in the last twenty years, using 
balloon-borne experiments (Leventhal et a!. 1977), and detectors on-board satellites such 
as HEAO-3 (Mahoney et al. 1982), SMM (Leising et al. 1988), and CGRO (Leising 1993; 
Iyudin et al. 1995), from which upper limits on the ejected 22 Na have been derived. In 
particular, the observations performed with the COMPTEL experiment on-board CGRO of 
five recent Ne-type novae (Nova Her 1991, Nova Sgr 1991, Nova Set 1991, Nova Pup 1991 
and Nova Cyg 1992. Iyudin et al. 1995), as well as observations of standard CO novae, 
have led to an upper limit of 3.7 x 10 -8 M for the 22 Na mass ejected by any nova in the 
Galactic disk. This restrictive limit has posed some constraints on pre-existing theoretical 
models of classical nova explosions. 

26 Al is another unstable nucleus, with a lifetime of r = 1.04 x 10 6 years, that decays 
from ground state to the first excited level of 26 Mg , which in turn de-excites to its ground 
state by emitting a gamma-ray photon of 1.809 MeV. This characteristic gamma-ray 
signature, first detected in the Galatic Center by the HEAO-3 satellite (Mahoney et al. 
1982, 1984), has been confirmed by other space missions like SMM ( [Share et al. 1983 ) 



and by several balloon-borne experiments. The most recent measurements made with 
COMPTEL have provided a map of the 1.809 MeV emission in the Galaxy (Diehl et al. 
1995, 1997; Prantzos & Diehl 1996). The inferred 1 - 3 M of Galactic 26 A1 are, according 
to the observed distribution, mainly attributed to young progenitors, such as massive AGB 
stars, type II supernovae and Wolf-Rayet stars. More recent analyses of the 1.809 MeV 
COMPTEL map reveal a correlation between this map and the COBE/DMR maps, tracing 
free-free emission, thus confirming that the main contributors to the Galactic 26 A1 are 
massive stars (Knodlseder 1997). However, a potential contribution from novae or low-mass 
AGB stars cannot be ruled out (see Jose, Hernanz & Coc 1997, for a recent analysis of 26 A1 
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production in classical novae). 

First estimates of the 22 Na and 26 A1 production in novae were performed by different 
groups, using simplified one-zone models with representative temperature and density 
profiles. Hillebrandt & Thielemann (1982) and Wiescher et al. (1986) suggested already 
that classical novae might produce significant amounts of 26 A1, not enough to represent 
major Galactic sources, but relevant to account for the observed isotopic anomalies found 
in some meteorites. New parametrized calculations by Weiss & Truran (1990) and Nofar, 
Shaviv & Starrfield (1991), revealed that nova envelopes previously enhanced in heavy 
elements (from Ne to Mg) produce large amounts of 22 Na and 26 A1. Since this metal 
enrichment is expected to result from dredge-up of core material, Weiss & Truran suggested 
that massive ONeMg white dwarfs (the ones attaining the highest peak temperatures 
and, therefore, the most efficient dregde-up) are likely to provide the largest abundances 
of both 22 Na and 26 A1 in the ejecta. Politano et al. (1995) revisited this scenario using 
hydrodynamic computations. They reported on a strong anticorrelation between 22 Na 
and 26 A1 production: novae that produce the largest amounts of 22 Na (i.e., massive white 
dwarfs) are not the same as those accounting for the largest yields of 26 A1 (i.e., low- mass 
white dwarfs). The obtained 22 Na yields range between 5 x 10~ 5 and 5 x 10~ 3 , by mass. 
Assuming that the whole envelope (~ 10~ 4 — 10~ 5 M Q ) is ejected during the outburst, 
they concluded that nearby ONeMg novae with M W d > 1.25 M should produce detectable 
22 Na 7-rays for CGRO, a prediction not confirmed so far (Iyudin et al. 1995). Their results 
showed also a significant production of 26 A1 (i.e., (19.6 — 7.5) x 10~ 3 , by mass, corresponding 
to ONeMg white dwarfs with masses between 1.0 — 1.35 M , respectively), which could 
account for a major fraction of the Galactic 26 A1. 

Recent hydrodynamic computations of ONe novae (Jose, Hernanz & Coc 1997; Jose & 
Hernanz 1997, 1998) using updated initial compositions and nuclear reaction rates, have 
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led to a significant reduction of both 26 A1 and 22 Na ejected during nova outbursts, fn 
particular, a mean mass fraction of 1 x 10~ 4 of 22 Na is found in the 1.25 M ONe Model 
(with M e j ec ( 22 Na) = 1.3 x 10~ 9 M ), whereas a maximum value of 6 x 10~ 4 results from 
the 1.35 M ONe Model (M eiec ( 22 Na)=2.6 x 10~ 9 M ). The corresponding peak fluxes 
in the 1.275 MeV 22 Na line, below 10~ 5 photons s _1 cm -2 for novae at 1 kpc, turn out to 
be too low to be detected with OSSE or COMPTEL but represent potential targets for 
the nearby future INTEGRAL mission (Hernanz et al. 1997; Gomez-Gomar et al. 1998). 
Concerning 26 A1, yields ranging from 2 x 10~ 3 to 2 x 10~ 4 by mass have been obtained in a 
series of ONe nova models with masses between 1.15 — 1.35 M . Contribution of novae to 
the Galactic 26 A1 is limited to 0.1 — 0.4 M (Jose, Hernanz & Coc 1997). Nevertheless, a 
larger contribution cannot be ruled out if the (uncertain) lower limit for ONe white dwarfs 
is reduced down to 1.0 M (Jose & Hernanz 1998). 

Other hydrodynamic computations performed by Starrfield et al. (1997, 1998), using 
also updated nuclear reaction rates and opacities, have modified their previous estimates 
(Politano et al. 1995). The expected abundance of 22 Na in the ejecta has risen up to 
(2 — 3) x 10~ 3 , by mass, when 1.25 M ONeMg white dwarfs are adopted, high enough to 
be detected by CGRO, provided that all the accreted envelope (3 x 10~ 5 M ) is ejected. 
On the other hand, the improved input physics translates into a factor of 10 reduction on 
the synthesis of 26 A1, in better agreement with the analysis of the 1.809 MeV emission map 
provided by COMPTEL, and also with the results previously reported by Jose, Hernanz & 
Coc (1997). 

Whereas the agreement between the different groups on the expected contribution 
of classical novae to the Galactic 26 A1 has significantly increased, there remains some 
discrepancy concerning the amount of 22 Na present in the ejecta. Since the synthesis of 
both 22 Na and 26 A1 is very dependent on the adopted nuclear reaction rates, the large 
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uncertainties present in some key reactions of both NeNa and MgAl cycles (Coc et al. 1995; 
Prantzos & Diehl 1996), may significantly modify the expected yields. In particular, the 
study of the influence of a given reaction rate on 22 Na production is not trivial due to the 
two possible modes of formation and their interplay with convection (Coc et al. 1995). 
Moreover, during a nova outburst, thermodynamic conditions change on a short timescale 
so that nuclear reactions are never close to equilibrium. Due to obvious experimental 
difficulties, reaction rates involving short-lived radioactive nuclei of the NeNa-MgAl group 
are in general poorly known and the associated uncertainties may reach many orders of 
magnitude. It is thus important to know how the yields of important isotopes like 22 Na 
and 26 A1 are affected by those uncertainties. In this paper, series of hydrodynamic nova 
models have been computed assuming upper, recommended and lower estimates of the 
reaction rates, from which limits on the production of both 22 Na and 26 A1 are derived. The 
present analysis is focused on capture rates on radioactive nuclei, and for a temperature 
domain in the range T§ ~ 0.5 - 3.5. We refer to the NACRE compilation (Angulo et al. 
1998) for a more general discussion concerning capture rates on stable isotopes, and a 
wider temperature range. For convenience, we use the terms NeNa cycle and MgAl cycle 
to denote the reactions involved in 22 Na and 26 A1 formation. However, because of the high 
23 Na(p,7) and 27 Al(p,7) rates, they cannot be considered as genuine cycles. 

In Section 2, we outline some details of the method of computation and input physics. 
A detailed analysis of the synthesis of 22 Na and 26 A1 in classical novae, together with 
the study of the role played by specific nuclear reactions of the NeNa-MgAl cycles, is 
given in Sections 3 & 4. Constraints on the production of 22 Na and 26 A1 assuming lower, 
recommended and upper rates for some key reactions are derived in Section 5. The most 
relevant conclusions of this paper are summarized in Section 6. A detailed Appendix 
focused on the nuclear physics aspects of the reaction rates within the NeNa and MgAl 
cycles follows. 



2. Model and Input physics 

Evolutionary sequences of nova outbursts have been calculated by means of an updated 
version of the code SHIVA (see Jose 1996; Jose & Hernanz 1998), a one- dimensional, 
implicit, hydrodynamical code in lagrangian formulation, that follows the course of the 
outburst from the onset of accretion up to the expansion and ejection stages. The code 
solves the standard set of differential equations for hydrodynamical evolution: conservation 
of mass, momentum and energy, energy transport by radiation and convection, plus the 
definition of the lagrangian velocity, including a time-dependent formalism for convective 
transport whenever the characteristic convective timescale becomes larger than the 
integration time step (Wood 1974). Partial mixing between adjacent convective shells 
is treated by means of a diffusion equation (see Prialnik, Shara & Shaviv 1979, for the 
formalism). The code is linked to a reaction network, which follows the detailed evolution 
of 100 nuclei, ranging from 1 H to 40 Ca, through 370 nuclear reactions, with updated rates, 
and screening factors from Graboske et al. (1973) and DeWitt, Graboske & Cooper (1973). 
As suggested by Politano et al. (1995), the matter transferred from the companion is 
assumed to be solar-like, and is mixed in a given fraction with the outermost shells of the 
underlying core by means of an unknown mechanism (either shear mixing, diffusion or a 
convective multidimensional process). The composition of the underlying core has been 
taken from recent detailed evolutionary models in the case of ONe white dwarfs, which are 
the main contributors to 22 Na and 26 A1 synthesis. These stars are made basically of 16 
and 20 Ne ( Pominguez, Tornambe fc Isern 1993J ; [Ritossa, Garcia-Berro fc Iben 1996| ), with 
smaller traces of 23 Na, 24 ' 25 Mg, 27 A1 and other species. This issue plays a crucial role in 
the resulting nucleosynthesis, and should be taken into account in order to compare results 
obtained by different groups. In particular, the ONeMg models computed by Starrfield et 
al. (1998) have an initial composition of the white dwarf core based on old nucleosynthesis 
calculations of C-burning from Arnett & Truran (1969), which is richer in 24 Mg and 20 Ne 
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than the one adopted in this paper (see Table 1). 

3. Synthesis of 22 Na in Classical Novae 

The high temperatures attained during ONe nova outbursts (with peak values within 
2 — 3.5 x 10 8 K), allow a noticeable nuclear activity in the NeNa and MgAl cycles, which 
results on a significant production of species of astrophysical interest, such as 22 Na and 
26 Al. In this Section, we will describe the main mechanisms of 22 Na synthesis, through 
a detailed analysis of a 1.25 M , ONe white dwarf, which accretes solar- like matter at a 
rate M = 2 x 10~ 10 M Q .yr _1 , assuming a 50% degree of mixing with the ONe core (Model 
ONe5, in Jose & Hernanz 1998). Snapshots of the evolution of several isotopes relevant to 
22 Na synthesis (i.e., 20,21,22 Ne, 21,22,23 Na and 22,23 Mg) are shown in Figure 1. 



3.1. Main nuclear reactions involved in 22 Na production 

At the onset of accretion, the evolution of 22 Na is mainly dominated by the chain 
of nuclear reactions 20 Ne(p, 7) 21 Na(,3 + ) 21 Ne(p, 7) 22 Na( / 9+) 22 Ne(p, 7) 23 Na (i.e., the cold 
mode of the NeNa cycle. See Fig. 2). When the temperature at the burning shell reaches 
Tig = 5xl0 7 K (Fig. 1, panel 1), the main nuclear reaction of the NeNa cycle is 21 Ne(p, 
7) 22 Na , which significantly reduces the amount of 21 Ne . In fact, when Tb s = 7xl0 7 K 
(panel 2), the amount of 21 Ne is already too small to maintain the main mechanism for 
22 Na synthesis. Therefore, 22 Na will begin to decrease near the burning shell due to proton 
captures, following the rise of temperature toward peak value (panels 3 & 4). 

At Tb s — 1x10 s K (panel 3), the amount of 21 Ne has already decreased below 10~ 6 
by mass, except at the outer envelope, where some 21 Ne is synthesized from the /9 + -decay 
of 21 Na (previously built up as 20 Ne(p, 7) 21 Na ). 21 Na and 23 Mg increase due to 20 Ne(p, 
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7) 21 Na and 22 Na(p, 7) 23 Mg , respectively. Destruction of 22 Na through (p, 7) reactions goes 
on but, due to convection, 22 Na shows a nearly flat profile throughout the envelope. 

When Tb s = 2xl0 8 K (panel 4), there is a dramatic decline in 23 Na (due to (p/y) 
and (p,a) reactions) and in the 22 Na abundance (by 22 Na(p, 7) 23 Mg ). 22 Ne has slightly 
decreased. Also noticeable is the increase of 21 Na (since 20 Ne(p, 7) 21 Na dominates 
destruction from both 21 Na(/3 + ) 21 Ne and 21 Na(p, 7) 22 Mg ), which plays a crucial role in the 
synthesis of 22 Na at the late stages of the outburst. Both 22 ' 23 Mg increase due to proton 
captures on 21,22 Na respectively. 

33 seconds later, the temperature at the burning shell attains its peak value, 
Tbs,max — 2.44 x 10 8 K (panel 5). The amount of 21 Na is maintained by a quasi-equilibrium 
between 21 Na(p, 7) 22 Mg and 20 Ne(p, 7) 21 Na . The mean amount of 22 Na increases due 
to 22 Mg(/3 + ) 22 Na , previously transported by convection to the outer, cooler layers of the 
envelope. The amount of 23 Na decreases due to both (p,7) and (p,a) reactions, which 
dominate 22 Ne(p, 7) 23 Na as well as 23 Mg(/9 + ) 23 Na . Since the peak temperature achieved 
in the burning shell is not extremely high, 20 Ne remains nearly unchanged (the typical 
temperature for 20 Ne burning exceeds ~4x 10 8 K). With respect to the other neon isotopes, 
21 Ne increases (mainly due to the /9 + -decay of 21 Na at the outer shells), whereas 22 Ne is 
destroyed by 22 Ne(p, 7) 23 Na . The two magnesium isotopes 22,23 Mg increase as a result of 
proton captures onto 21,22 Na. 

Shortly after, due to the sudden release of energy from the short-lived /9 + -unstable 
nuclei 13 N, 14,15 0, and 17 F, the envelope begins to expand. The role played by (p,7) and 
(p,a) reactions is therefore reduced following the drop in temperature, whereas /5 + -decays 
progressively dominate the evolution. The abundances of 21 Na and 22,23 Mg decrease as a 
result of such /3 + -decays, which in turn increase the amount of 21 Ne and 22 ' 23 Na (panels 6 
&7). 
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At the final stages of the outburst (panel 8), as the envelope expands and cools down, 
most of the remaining nuclear activity in the NeNa cycle is due to /3 + -decays, such as 
21 Na(/3+) 21 Ne , 22 Mg(/3+) 22 Na or 23 Mg(/3+) 23 Na . The resulting mean abundance of 22 Na 
in the ejected shells of this 1.25 M ONe Model is X( 22 Na ) = 9.6xl0~ 5 by mass, which 
corresponds to 1.3xl0 -9 M Q of 22 Na ejected into the interstellar medium. Other species 
of the NeNa cycle present in the ejecta are 20 Ne (X( 20 Ne ) = 0.18, a mass fraction slightly 
higher than its initial value because of the operation of the 23 Na(p,a) 20 Ne reaction. See 
Table 1), 21 ' 22 Ne (X( 21 Ne ) = 3.5 xKT 5 , and X( 22 Ne ) = 1.0 xKT 3 ), and 23 Na (with 
X( 23 Na) = 1.4xl0~ 3 ). 

3.2. Effect of the reaction rates on the synthesis of 22 Na 

From the abovementioned nuclear physics viewpoint, the synthesis of 22 Na is mainly 
controlled by four reactions: 20 Ne(p, 7) 21 Na , 21 Na(p, 7) 22 Mg , 21 Ne(p, 7) 22 Na and 22 Na(p, 
7) 23 Mg (Fig. 2), for which recent updates to the reaction rates are available. Despite 
energy production during nova outbursts is not very dependent on the specific prescriptions 
adopted for such rates (see Section 5.1), they play a crucial role on the accompanying 
nucleosynthesis, since changes in the reaction paths are expected. Several test models have 
been computed to analyse the role played by different nuclear reactions on the synthesis 
of 22 Na (as well as 26 A1). Main results concerning 22 Na and 26 A1 production, compared 
with the abundances found in Jose & Hernanz (1998) with previous prescriptions, are 
summarized in Table 2. 

20 Ne(p, 7) 21 Na is the slowest proton capture reaction on any stable neon and sodium 
isotope. According to the adopted composition for the ONe core (Ritossa, Garcia-Berro 
& Iben 1996), 20 Ne is the most abundant Ne-Na isotope. Therefore, 20 Ne(p, 7) 21 Na limits 
22 Na production in classical novae. The adopted 20 Ne(p, 7) 21 Na rate is the Caughlan & 
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Fowler (1988) one, which is based on Rolfs et al. (1975) data. In the considered domain of 
temperature, it is well known. The estimated uncertainty is of a factor of ~3 (defined as 
the ratio between high and low estimates within the domain of temperature considered), 
according to the new compilation of nuclear reaction rates (Angulo et al. 1998), but more 
likely a factor of ~1.5 only, as it is derived from short range standard extrapolations of 
experimental data (Rolfs et al. 1975). Therefore, nuclear uncertainties associated with 
20 Ne(p, 7) 21 Na should play no significant role on 22 Na production. 

Concerning the 21 Ne(p, 7) 22 Na rate, no significant effect of the reduction of the 
contribution of the E x = 6.384 MeV level (Gorres et al. 1983) with respect to Caughlan & 
Fowler (1988) has been found on 22 Na production (tested with a 1.25 M Q ONe white dwarf 
model). According to Angulo et al. (1998), other nuclear uncertainties scarcely affect the 
21 Ne(p, 7) 22 Na rate at moderate temperatures. On the contrary, much more uncertain is 
the 22 Ne(p, 7) 23 Na rate, poorly known in the range of temperatures of interest for nova 
outbursts (Angulo et al. 1998). Nevertheless, due to the low initial 22 Ne abundance with 
respect to 23 Na, and also to the negligible 22 Na(/9 + ) 22 Ne decay (r = 3.75 yr) during the 
TNR, the nuclear uncertainties affecting the 22 Ne(p, 7) 23 Na rate are not relevant for 22 Na 
production in novae. 

Recent experimental investigations of the 22 Na(p,7) 23 Mg reaction (Seuthe et al. 
1990; Schmidt et al. 1995; Stegmuller et al. 1996) have provided a firmer basis for the 
determination of its rate. Since the new rate is lower than the old Caughlan & Fowler's 
(1988) analytic fit (i.e., one order of magnitude for T§ > 1), destruction of 22 Na by means 
of (p,7) reactions is reduced. In a test model, consisting of a 1.25 M ONe white dwarf 
(Model ONe5, Jose & Hernanz 1998), the mean abundance of 22 Na in the ejecta increases 
by a factor of ~ 3 when the new 22 Na(p, 7) 23 Mg rate (Stegmuller et al. 1996), instead 
of the CF88 one, is adopted (see Table 2). It is worth noticing that significant nuclear 
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uncertainties affect this rate (a factor ranging from 3 to 6 for T 8 > 1), which turn out to be 
crucial in order to derive ranges of 22 Na production during nova outbursts. According to 



Stegmuller et al. (1996)| , this uncertainty is mainly due to the possible effect of a resonance 
at E p = 225 keV. Below T§=1, the uncertainty reaches three orders of magnitude but the 
rate remains small enough to prevent destruction of 22 Na. 

Three test models have been computed to analyze the role played by 21 Na(p, 7) 22 Mg , 
in view of the nuclear uncertainties present in this rate, in particular the estimated strength 
of the first E x = 5.714 MeV level above the proton threshold (see Appendix, section A.l). 
Calculations assume ONe white dwarfs of masses 1.15, 1.25 and 1.35 M Q , and the same 
input physics than Models ONe3, ONe5 and ONe6 described in Jose & Hernanz (1998), but 
reducing the 21 Na(p, 7) 22 Mg rate given by CF88 by a factor of 100 (similar, for novae, to the 
lower rate given in the Appendix). It results in a significant increase in the 22 Na abundances 
present in the ejecta (a factor of ~ 2 to 3, in the 1.15 and 1.25 M Models, and a factor of ~ 
1.2, in the 1.35 M Model. Table 2), as compared with the values found with the standard 
CF88 rate. This effect can be interpreted as follows: when the 21 Na(p,7) 22 Mg(/3 + ) 22 Na rate 
is reduced by a factor of 100, the alternative path, 21 Na(/9 + ) 21 Ne(p,7) 22 Na, is favored. In 
this case, 22 Na production is delayed to a time when the envelope is already expanding and 
cooling down (contrary to the case when the higher 21 Na(p, 7) 22 Mg rate is adopted). As 
a result, a major fraction of 22 Na survives. This, in turn, explains the lower effect found 
in the 1.35 M Model, caused by the higher temperatures achieved in the envelope (with 
Tpeak = 3.2 x 10 8 K) which remain high enough at the time when 22 Na is synthesized. One 
should note that this effect (increase in the 22 Na yield when the 21 Na(p,7) 22 Mg rate is 
reduced) was not foreseen and stresses the importance of full hydrodynamical calculations. 

Other reactions that may be involved in the synthesis of 22 Na (and 23 Na) are 
22 Mg(p,7) 23 Al, 23 Al(p,7) 24 Si, or 23 Mg(p,7) 24 Al. However, no significant effect on the 22 Na 
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production is found when using updated rates for such reactions (see Appendix and Table 
2), because of the limited nuclear flow they conduct for nova conditions. 

According to this analysis, full evolutionary sequences of nova outbursts, from the 
onset of accretion up to the ejection stage, have been performed taking into account the 
low, recommended and high estimates to the 21 Na(p, 7) 22 Mg , 22 Na(p, 7) 23 Mg and 21 Ne(p, 
7) 22 Na rates (as well as the new 23 Mg(p,7) 24 Al rate derived by Kubono, Kajino & Kato 
1995), in order to derive a range of 22 Na yields, resulting purely from nuclear physics 
uncertainties. Results are summarized in Section 5. 

4. Synthesis of 26 A1 in Classical Novae 

In this Section, we will analyze the nuclear paths leading to 26 A1 synthesis. As for 22 Na 
production, we will describe the course of the 1.25 M , ONe nova model (Jose & Hernanz 
1998, Model ONe5). Snapshots of the evolution of several isotopes relevant to 26 A1 synthesis 
(i.e., 24 > 25 > 26 Mg, 25 ' 27 A1, 26 ' 27 Si and the ground and isomeric states for 26 A1, herefater 26 A1 9 
and 26 Al m ) are shown in Figure 3. 

4.1. Main nuclear reactions involved in 26 A1 production 

Nucleosynthesis of 26 A1 is complicated by the presence of a short lived {jxn = 6.3 s) 
spin isomer. At low temperatures (T 8 ^ 4), the 26 A1 ground and isomeric states do not 
reach thermal equilibrium and must be treated as two separate isotopes (Ward & Fowler 
1980, see also Coc & Porquet 1998). 

The nuclear activity in the MgAl cycle at the early phases of the TNR, when the 
temperature at the burning shell is Tb s = 5 x 10 7 K, is dominated by 25 Mg(p,7), which 
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leads to both 26 A1 ground and isomeric states. A significant amount of 26 A1 9 is already 
synthesized at such temperatures (Figure 3, panel 1). Another aluminum isotope, 27 A1, is 
slightly enhanced with respect to its initial abundance by means of 26 Ar m ( / 5 + ) 26 Mg(p,7) 27 Al. 
A similar trend is found when T\, s = 7 x 10 7 K (panel 2): both 26 A1 9 and 27 A1 continue to 
rise (specially the first one, which increases by nearly a factor of 10 from the abundance 
shown in panel 1). 

At Tb s = 10 8 K (panel 3), the evolution in the MgAl cycle is dominated by 
24 Mg(p,7) 25 Al(/9 + ) 25 Mg, which in turn accounts for a noticeable production of the 
/3 + -unstable nuclei 25 A1, and also for the increase in the mean 26 A1 9 abundance (by means 
of 25 Mg(p,7) 26 Al 9,m ). A nearly flat profile of 26 A\ 9 results from convective mixing, which 
extends already throughout the whole envelope. 

A major change in the dominant nuclear path is found when Tj, s = 2 x 10 8 K (panel 
4). The temperature attained near the burning shell is high enough to allow (p/y) reactions 
to proceed efficiently In particular, the abundance of 24 Mg is reduced by a factor of 
~ 1000. The isomer 26 Al m exceeds already 10~ 3 by mass for most of the envelope. A 
significant fraction of 25 A1 is transformed through proton captures into 26 Si, which will 
increase the abundance of 26 Al m , and in turn that of 27 A1 in the late phases of the TNR. 
The final 27 A1/ 26 A1 9 ratio will reflect a competition between two different paths: 24 Mg (fed 
by 23 Na(p,7) 24 Mg) is transformed by proton captures into 25 A1, which can either decay 
into 25 Mg or capture another proton, leading to 26 Si. Only the first channel accounts for 
26 Al 9 synthesis, whereas 27 A1 can be produced by both paths, following 25 Mg(p,7) 26 Al fl,m (p, 
7) 27 Si( / 3+) 27 Al, or 26 Si(/3 + ) 26 Al m (/3+) 26 Mg(p, 7 ) 27 A1. Anyway, because of the large 
abundance of 25 Mg in the envelope at this stage, proton captures on the seed 25 Mg become 
the major source of 26 A1 9 . Some leakage from the MgAl cycle due to 26 Al s ' m (p,7) 27 Si is also 
obtained at this stage of the outburst. 
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When the burning shell attains its peak temperature (T{, Sjmox = 2.44 x 10 8 K. Panel 
5), most of the MgAl isotopes show a significant reduction near the burning shell because 
of (p,7) reactions. The dominant paths at this stage are 25 Mg(p,7) 26 Al 9 (p,7) 27 Si and 
27 Al(p,7) 28 Si, which also account for a significant leakage from the cycle (moreover, 
26 ' 27 Si(p,7) are much faster than the corresponding /3 + -decays). Following the course of 
the outburst (panels 6 to 8), as the envelope expands and cools down, proton capture 
reactions are progressively reduced. Therefore, the late time evolution is mainly dominated 
by /3+-decays, such as 26 Si(/3+) 26 Al m , 27 Si(/3+) 27 Al, 25 Al(/3+) 25 Mg and 26 Al m (/3+) 26 Mg, 
which in turn, increase the final amounts of 25 ' 26 Mg and 27 A1. The mean abundance of 26 A1 
in the ejecta of this 1.25 M , ONe white dwarf model is X( 26 A1) = 5.4 x 10~ 4 by mass, 
which translates into 7.6 x 10~ 9 M of 26 A1 ejected into the interstellar medium. Other 
isotopes of the MgAl group present in the ejecta are Al (X( 27 A1) = 2 x 10~ 3 , half the 
initial abundance. See Table 1), 25 Mg (2.4 x KT 3 ), 26 Mg (2.8 x KT 4 ), and 24 Mg (2 x 1(T 4 ). 

4.2. Effect of the reaction rates on the synthesis of 26 A1 

The previous analysis has revealed that several isotopes should be considered as 
potential seeds for 26 A1 synthesis: 24 Mg , 25 Mg and, to some extent, 23 Na and 22 Ne . In 
this case, the number of nuclear reactions involved in the synthesis of 26 A1 is rather large 
(see Fig. 2). Several test models of nova outbursts have been computed to analyze the role 
played by each one of the most relevant reactions. Main results are summarized in Table 2. 

The 23 Na(p,7) 24 Mg rate (in contrast with 23 Na(p,a) 20 Ne) has been strongly modified 
by the introduction of the upper limit for the E^. m =0.138 MeV resonance strength obtained 
by Gorres, Wiescher & Rolfs (1989). When taking this upper limit, the rate is significantly 
increased up to T 8 ~2 (Gorres et al. 1989). Its effect has been tested by recomputing a 1.15 
M , ONe white dwarf Model (i.e., Model ONe3. Jose & Hernanz 1998): an increase by 
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a factor of ~ 3 of the final 24 Mg yields results, which in turn leads to slightly larger mass 
fractions of both 25 Mg (factor 1.4) and 26 A1 (factor 1.3). Some 27 A1 is also overproduced 
(by a factor of 1.4) when this resonance is taken into account. Those effects are expected 
to be stronger for more massive white dwarfs. 

At the temperatures attained in nova outbursts, 24 Mg(p,7) 25 Al may become even 
faster than 12 C(p,7) 13 N. The 24 Mg(p,7) 25 Al rate is dominated by the presence of the E^ m 
= 0.214 MeV, J 7r = l/2 + resonance and suffers little uncertainty. 

25 Mg(p,7) 26 Al 9 is also of great importance since it is the only channel that leads to 26 A1 
in its ground state. One may expect that its potential uncertainties are directly reflected in 



the 26 A1 9 yields. The new rate is significantly smaller by a factor of ~7 ([Hiadis et al. (1996)|) 



than the |Caughlan fc Fowler 1988| one, below Tg ~ 1, but this has no consequences in the 
nova domain, lliadis et al. (19975]] reanalyzed available transfer reaction data concerning the 



lower lying resonances and concluded that the associated uncertainties are relatively small 
(i.e., a factor of 1.5-2). For T 8 > 1.5, available direct measurements lead also to very small 
uncertainties (a factor of 1.5, according to |Iliadis et al. (1996)|) and should not influence 



significatively 26 A1 production. Moreover, since 25 Mg(p,7) can lead either to 26 A1 9 or 26 Al m 
with a known branching ratio (Endt & Rolfs 1987, lliadis et al. 1996), the 25 Mg(p,7) 26 Al m 
rate is known with a similar precision as the 25 Mg(p,7) 26 Al 5 one. 

In nova nucleosynthesis, the most important resonance in 26 Al s (p,7) 27 Si is the E^ m = 
0.188 MeV, while uncertainties on the lower lying resonances have no influence (Coc et al. 
1995). In our calculations, we use the only single direct measurement available (Vogelaar 
1989) for this crucial resonance strength. Recently, the results of the proton transfer 
reaction that motivated this subsequent direct measurement have been published (Vogelaar 
et al. 1996) but the resonance strength deduced from this indirect measurement is strongly 
dependent on the adopted transferred angular momentum. Accordingly, we have checked 
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the influence of such resonance on the resulting 26 A1 yields by reducing its strength by a 
factor of 1/3, well within the range of values adopted by [Angulo et al. (1998)| . We obtain a 



significant increase by a factor of ~ 2 in the final amount of 26 A1 (for a 1.15 M ONe white 
dwarf Model. Table 2), which strongly stresses the need of additional direct measurements 
to confirm the results found by Vogelaar (1989). 

25 Al(p,7) 26 Si plays an important role on 26 A1 synthesis, since it leads to the formation 
of the short-lived isomer (through 26 Si(/3 + ) 26 Al m ), instead of the long-lived ground state. 
We have checked the role played by the 25 Al(p,7) 26 Si rate, adopting the upper and lower 
estimates (see Appendix) provided by Coc et al. (1995). Results are compared with those 
obtained with the 25 Al(p,7) 26 Si rate given by Wiescher et al. (1986). Whereas the final 
amounts of 26 A1 and 27 A1 are very slightly enhanced when the lower rate is adopted (i.e., 
case A. Coc et al. 1995), a significant reduction of 26 A1 by a factor of ~ 2 is obtained for 
the upper one (case C). In the first case, 25 Al(/5 + ) 25 Mg becomes faster than 25 Al(p,7) 26 Si 
and, therefore, the final amount of 25 Mg (and in turn 26 A1 ) increases (Table 2). 

Large uncertainties affect the 26 Al m (p,7) 27 Si rate. The existing Caughlan & Fowler's 
(1988) prescription for this rate results from a Hauser-Feschbach calculation, a statistical 
approach that is not reliable at low temperatures, where the contribution of isolated 
resonances dominates. Hence, the rate may differ by several orders of magnitude from the 
theoretical one at the temperatures achieved during nova outbursts. However, no noticeable 
effect on 26 A1 production is found when the 26 Al m (p,7) 27 Si rate, as given by CF88, is 
multiplied arbitrarily by a factor of 100, with the exception of a net reduction (factor ~ 2) 
on the final amount of 26 Mg . 

Nuclear uncertainties significantly affect the 26 Mg(p,7) 27 Al rate. Since the initial 
amount of 26 Mg is of the same order of magnitude than that of 27 A1, and it is also fed 
continuously by 26 Al m (/3 + ) 26 Mg (r = 9.15 s), these uncertainties cannot be ignored as for 
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22 Ne(p,7) 23 Na (see Section 3.2). They are specially important around Tg = 0.5, where 
the rate remains uncertain by a factor of ~ 100 (Champagne et al. 1990), due to the 
unknown strength of an hypothetical resonance at E c ^ = 90 keV. No influence on the final 
26 ' 27 A1 yields results from the inclusion of the 26 Mg(p,7) 27 Al recommended rate given by 
Champagne et al. (1990) (tested with a 1.15 M ONe white dwarf Model). 

With respect to CF88, the 27 Al(p,7) 28 Si rate has not changed: above T 8 =l, 
uncertainties are limited to a factor of ~2 (Angulo et al. 1998). On the contrary, new 
experimental data for the 27 Al(p,a) 24 Mg rate appeared shortly after the CF88 compilation: 
a direct measurement by [limmermann et al. (1988)| and a study of proton and alpha 
emission from 28 Si levels performed by |Uhampagne et al. (1988")| . As a result, the rate 
is strongly reduced with respect to the Caughlan & Fowler (1988) one by up to 4 orders 
of magnitude in the region of interest. As a consequence, the calculated 26 A1 yields were 
strongly reduced (Jose, Hernanz & Coc 1997). 

The new rates available for 23 Al(p, 7 ) 24 Si (Schatz et al. 1997), 26 Si(p,7) 27 P (Herndl 
et al. 1995), and 23 Mg(p, 7 ) 24 Al (Kubono, Kajino & Kato 1995 and Herndl et al. 1998) 
have no significant effect on the resulting 26 A1 and 27 A1 yields, as compared with the 
results obtained by Jose & Hernanz (1998) with 1.25 M ONe white dwarfs, using earlier 
prescriptions for these rates (i.e., van Wormer et al. 1994, Wiescher et al. 1986, and also 
Wiescher et al. 1986, respectively). 

5. Updated nuclear reaction rates, nuclear uncertainties and their effect on 

the synthesis of 22 Na and 26 A1 

Constraints on the synthesis of 22 Na and 26 A1 during nova outbursts have been derived 
by means of hydrodynamic calculations involving ONe white dwarfs of masses 1.15, 1.25 and 
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1.35 M , assuming lower, recommended and upper estimates for several key reactions of 
the NeNa-MgAl cycles. Table 3 lists the set of reaction rates adopted for the calculation of 
minimum to maximum 22 Na and 26 A1 production, compared with the rates adopted in Jose 
& Hernanz (1998). The specific choice adopted for these three sets is directly determined 
by the range of nuclear uncertainties accompanying the rates. The corresponding yields 
obtained with the different nuclear reaction networks, which are summarized in Tables 4 to 
6, allow to derive error bars to the synthesis of 22 Na and 26 A1. Since we have not included 
other sources of uncertainty (convection, modelization of the explosion, ...) these error bars 
are only of nuclear physics origin. 

5.1. Characteristics of the explosion 

The early stages of the outburst are mainly dominated by the CNO cycle (Jose & 
Hernanz 1998). Therefore, the update of the nuclear reaction rates of the NeNa-MgAl cycles 
has no influence on the characteristics of the accretion phase (for instance, the duration of 
the accretion phase, t acc , or the mass of the accreted envelope, AM env ). 

Differences in the time evolution appear when the temperature near the burning shell 
reaches T^ s ~ 10 8 K. At this stage, energy generation by nuclear reactions involves a 
relevant contribution from the NeNa-MgAl cycles, together with the hot and cold modes of 
the CNO cycle. The use of updated rates modifies several properties of the TNR, such as 
the time required for a temperature rise from 3 x 10 7 K up to 10 8 K (hereafter, t r i se ). For 
instance, a value of t r i se = 6.8 x 10 6 s was found for Model ONe5 (Jose & Hernanz 1998), 
whereas a shorter time, t r i se = 4.3 x 10 6 s, has been obtained in Model ONel25B, computed 
with updated (recommended) rates. 

The role played by the NeNa-MgAl cycles increases when the temperature at the 
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location of the burning shell exceeds 2 x 10 8 K. This has an important effect on the energy 
production as well as on the resulting peak temperature achieved during the TNR. For 
instance, peak values for the nuclear energy generation rate and temperature achieved in 
Model ONe5, e nuC)max = 2.1 x 10 16 ergg^s" 1 and T max = 2.44 x 10 8 K, translate into 
e n uc,max = 3.0 x 10 16 ergg^s" 1 and T max = 2.51 x 10 8 K, for Model ONel25B. A similar 
trend was also pointed out by Starrfield et al. (1998), using updated nuclear reaction rates 
with respect to Caughlan & Fowler's (1988) ones. 

Such differences have also some influence on the last phases of the evolution, when 
the envelope expands, cools down and eventually a fraction of the formerly accreted shells 
is ejected. In general, we find that models computed with updated rates lead to slightly 
larger ejected masses with larger mean kinetic energies (see Tables 4-6) as a result of the 
larger e nuc>max attained. It is worth noticing that these differences in the properties of the 
TNRs will be reflected in the accompanying nucleosynthesis. The reason is twofold: first, 
the update of the reaction rates implies a certain modification of the nuclear paths. And 
second, the differences found in the time evolution as well as on the peak temperatures 
achieved modify the role played by (p,7) and (p,a) reactions. 

5.2. Nucleosynthesis in the NeNa-MgAl cycles 

Due to the higher peak temperatures achieved in the models computed with the 
updated network (recommended rates. See Tables 4-6), 20 Ne is slightly reduced in the 
ejecta with respect to previous results obtained with older rates (Jose & Hernanz 1998). 
Differences are rather small (i.e., a mean mass fraction of 0.17 by mass, instead of 0.18, 
results for the 1.15 & 1.25 M models), but they turn out to be important to determine the 
synthesis of less abundant species within the NeNa cycle. For instance, the abundance of 
21 Ne increases by a factor between ~ 3 — 5 (the higher temperatures as well as the lower rate 
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adopted for 21 Na(p,7) 22 Mg favor the chain 20 Ne(p,7) 21 Na(/3 + ) 21 Ne), however no significant 
change is found for 22 Ne (which is essentially reduced by (p,7) reactions from its initial 
amount). Worth noticing is also the net increase in the final amount of 22 Na in the ejecta: 
a factor of ~ 4 — 5 for the 1.15 & 1.25 M Q models (resulting from the lower rates adopted 
for 21,22 Na(p,7) 22,23 Mg). This may have important consequences for gamma-ray astronomy, 
since it translates into a change by a factor of ~ 2 in the maximum expected distance at 
which the 1275 keV 22 Na line emitted by an exploding ONe novae would be eventually 
detected. Other isotopes, such as 23 Na and 25 Mg, are also overproduced (except for the 
1.35 M model), but less efficiently As a result, the isotopic ratios 23 Na/ 22 Na decrease by 
nearly a factor of ~ 2 — 3. The abundances of the other magnesium isotopes, 24 ' 26 Mg, are 
slightly enhanced (in some cases even by a factor of ~ 2 — 4). Another interesting feature is 
that 26 A1 remains essentially unaffected by the update of the network when the recomended 
rates are used. Therefore, the conclusions relative to the small contribution of classical nova 
outbursts to the synthesis of the Galactic 26 A1 (see Jose, Hernanz & Coc 1997) still hold. 
Since the amount of 27 A1 (and that of 28 Si) in the ejecta remains essentially unchanged, no 
variation on the isotopic ratio 26 A1/ 27 A1 is found. 

Another relevant outcome from the nucleosynthetic viewpoint is the dispersion in the 
mean ejected abundances obtained when the uncertainties associated with the nuclear 
reaction rates are taken into account. In the following, we will analyse the impact of such 
uncertainties in the resulting mean abundances in the ejecta. For that purpose, we will 
compare the yields obtained when the two extreme networks listed in Table 3 (i.e., A and 
C, leading to maximum and minimum 22 Na- 26 Al production, respectively) are adopted. 
We define the dispersion factor, F, as the ratio between the mean ejected abundances of a 
given nuclear species obtained when network A and C are adopted. All models computed 
(i.e., involving 1.15, 1.25 or 1.35 M white dwarfs) show a very small dispersion in the 
20 Ne yields, which results essentially from the different peak temperatures attained at the 
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envelope rather than from uncertainties in the nuclear reactions (such as 20 Ne(p,7) 21 Na). 
No relevant dispersion is found either for 22 Ne, 23 Na, 27 A1, and 28 Si. On the contrary, 
the 21 Ne yields show a rather wide dispersion: a factor F ~ 25 in the 1.35 M Q models. 
Concerning 22 Na, differences between F ~ 2 — 4 are found. Similar degrees of dispersion 
are also obtained for the magnesium isotopes, with F ranging from 3 to 5 in the case of 
25 Mg and from 2 to 3 for 24 Mg. Dispersion factors between 0.3 and 0.7 result for 26 Mg. The 
variation in the ejected amounts of 26 A1 is particularly worth noticing, with a dispersion 
factor ranging from F ~ 4 to 7 (not considering the possible uncertainty associated with 
the 188 keV resonance in 26 Al 9 (p,7) 27 Si). Also interesting to notice is the fact that the 
26 A1/ 27 A1 ratio remains in the range 0.5-0.1 (see Tables 4-6). 

In view of the abovementioned analysis, the nuclear uncertainties accompanying the 
reaction rates within the NeNa-MgAl cycles introduce a relatively wide dispersion in the 
yields resulting from classical nova outbursts. Therefore, predictions of 22 Na and 26 A1 yields 
would benefit from new nuclear physics experiments aimed at reducing the uncertainties 
associated with some key reactions of the NeNa-MgAl cycles, in particular 21 Na(p, 7) 22 Mg 
and 25 Al(p,7) 26 Si whose rates are uncertain by several orders of magnitude but also 
23 Na(p,7) 24 Mg and 22 Na(p, 7) 23 Mg . In addition, a verification of the yet unpublished 
values corresponding to the 0.188 MeV resonance of 26 Al 5 (p,7) 27 Si , as measured by 
Vogelaar (1989)| , is also recommended, due to its crucial role on the synthesis of 26 A1. These 



experiments should be focused on the measurement of a few key parameters (level energies, 
spectroscopic factors, ...). Many of these measurements involve short lived radioactive 
species and should benefit from the current development of radioactive ion beam facilities. 

On the contrary, our calculations have shown that the remaining uncertainties on 
22 Mg(p, 7 ) 23 Al, 23 Al(p, 7 ) 24 Si, 23 Mg(p, 7 ) 24 Al, 24 Al(p, 7 ) 25 Si, 26m Al(p, 7 ) 27 Si, 26 Si(p, 7 ) 27 P, 
26 Mg(p,7) 27 Al, 27 P(p,7) 28 S or 27 Si(p,7) 28 P have no effect on 22 Na and 26 A1 production. 
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6. Conclusions 

We have computed a series of hydrodynamic models of nova outbursts, from the onset 
of accretion up to the ejection stage, for a range of ONe white dwarfs with masses between 
1.15 and 1.35 M , with the aim of analyzing in detail the main nuclear paths leading to the 
synthesis of 22 Na and 26 A1. The role played by several key reactions within the NeNa-MgAl 
cycles has been tested by a series of models, which have been compared with previous 
results obtained by Jose & Hernanz (1998) with different prescriptions for the reaction 
rates. Limits on the production of both 22 Na and 26 A1 have been derived from a series of 
computations assuming upper, recommended or lower estimates to the reaction rates. The 
most relevant conclusions extracted from this work can be summarized as follows: 

1. The update of the nuclear reaction network results in a net increase in the final 
amount of 22 Na ejected into the interstellar medium during classical nova outburts. This 
translates into an increase (by a factor of ~ 2) of the expected maximum distance at which 
an exploding ONe nova would be eventually detected through its emission at 1275 keV 
( 22 Na decay 7-ray line). 

2. The final amount of 26 A1 remains essentially unaffected by the update of the 
network, confirming that classical novae scarcely contribute to the Galactic 26 A1, as pointed 
out earlier by Jose, Hernanz & Coc (1997). 

3. Large nuclear uncertainties affect some key reactions of the NeNa-MgAl cycles, 
with a significant effect on the production of both 22 Na and 26 A1. When a combination of 
reaction rates, leading to maximum or minimum 22 Na- 26 Al synthesis, is adopted, a large 
dispersion in the final abundances is found. We stress that the derived ranges for 22 Na- 26 Al 
production can be interpreted as error bars on sodium & aluminum production posed by 
nuclear physics uncertainties. 
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4. In order to reduce the impact of the nuclear uncertainties in the production of 
22 Na and 26 A1, we point out some nuclear reactions that deserve new experiments, in 
particular, 21 Na(p, 7) 22 Mg , 22 Na(p, 7) 23 Mg and 25 Al(p,7) 26 Si . A confirmation of the values 
corresponding to the 0.188 MeV resonance of 26 Al 5 (p,7) 27 Si , as measured by |Vogelaar 
|(1989)| , would also be of great importance. 
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A. Appendix: Nuclear reactions 

We discuss here the reaction rates corresponding to proton capture reactions on 
short-lived nuclei. Other reactions within the NeNa and MgAl cycles are analyzed in the 
text (see Sections 3.2 and 4.2) and more extensively in the forthcoming Compilation of 
Charged-P article Induced Thermonuclear Reaction Rates (Angulo et al. 1998). First, we 
introduce some standard nuclear physics quantities and notations to be used hereafter. 
We also briefly summarize several indirect methods that are generally used to extract or 
estimate unavailable direct data, for the calculation of thermonuclear rates. 

For an A(x,y)B reaction, resonance strengths (at E = E r ) are given by wy, where 
uj = (2J r + 1) ((2J X + 1)(2Ja + 1))~ is the statistical spin factor and 7 = T x T y /T 
is the width ratio. The total width (r(E')) corresponds to the sum of partial widths 
(Ti(E)) : T = J2i= x ,y,... Ft- The spectroscopic factor, C 2 S, relates single particle widths 
(T x p -) to actual ones: T x = C 2 S T x p \ These single particle widths are calculated 
by solving the Schrodinger equation for a particle in a realistic nuclear potential. An 
approximation for the upper limit of single particle widths is given by the Wigner limit: 
rjf = 3h 2 (fis 2 )' 1 Pe x (E = E r ,s), where s is the channel radius, \x the reduced mass, 
Pe x (E,s) the penetration factor and £ x the orbital angular momentum transferred by 
particle x. We introduce the reduced width, 9 2 , defined by T x = d 2 x Y^ , as it is sometimes 
used to derive estimates of reaction rates (see Iliadis 1997 for a dicussion of the relationship 
between Q 2 and C 2 S). 

Experimental spectroscopic factors can be extracted from transfer reactions (e.g., 
the proton transfer reaction X( 3 He,d)Y will lead to one proton spectroscopic factors 
used to calculate proton widths entering into the calculation of the X(p/y)Y rate). 
Spectroscopic factors for conjugate reactions (obtained by the p v^ n exchange) are 
assumed to be approximately equal. For instance, neutron spectroscopic factors for 22 Ne 
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can be extracted from the one neutron 21 Ne(d,p) 22 Ne transfer reaction. According to the 
above approximation, these neutron spectroscopic factors can be used to determine proton 
widths for the calculation of the 21 Na(p,7) 22 Mg rate for instance. Radiative widths (r 7 ) 
corresponding to transitions of the magnetic (MA) or electric type (EX) of order A are 
often expressed in Weiskopf's units (W.u.) to remove the effect of their i? 2A+1 dependence 
and strong variation with A (e.g. Firestone et al. 1996). Compilations of T 7 , expressed 
in these units, are available (Endt 1979) and may be used to estimate unknown radiative 
widths. They can also be obtained from the conjugate nuclei after correction for the 
i? 2A+1 dependence. If particle emission is energetically forbbiden in the conjugate level, the 
radiative width equals the total width. If the lifetime of this level is known, the radiative 
width is easily obtained. 

A.l. The 21 Na(p,7) 22 Mg reaction 



Estimates of the 21 Na(p,7) 22 Mg reaction rate QUaughlan fc Fowler 1988 ) have been 



provided by [Wiescher et al. (1986)| and [Wiescher fc Langanke (1986)| , considering the 



first three levels (Endt 1990) above the proton threshold. With the exception of the first 
level, the proton widths are much larger than the gamma widths, so that wy p» ooT t Only 
the first two levels (E x =5.714 and 5.837 MeV) contribute to the rate in the temperature 
domain considered. The E x = 5.837 MeV level is assumed to be the conjugate of the E x 
= 5.910 one in 22 Ne (Endt 1990). If this assignment is correct, the corresponding gamma 
width can be reliably deduced (Wiescher & Langanke 1986). On the contrary, the strength 
of the first E x — 5.714 MeV, J 7T =2 + level suffers from a significant uncertainty. For this 
level, the total width is known experimentally to be r=16.5±4.4 meV. To calculate the 
corresponding resonance strength, |Wiescher fc Langanke (1986)| and |Wiescher et al. (1986) 



estimated the proton width, r„, and deduced the radiative width, T 7 , from the relation 
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r = T p + r~. To estimate the proton width, two hypotheses were considered: |Wicschcr 



et al. (1986)| assumed £ p = and ffi = 0.01, but later [Wiescher fc Langanke (1986)| took 



£ p = 2 and ft; = 0.5. Nevertheless, because of the T p + T 7 = 16.5±4.4 meV constraint, the 
corresponding width ratios, 7, are not much different: 3.4 or 3.8 meV, very close to the 
maximum value (r/4), obtained when T 7 = T p = T/2. There are two reasons to think that, 
on the contrary, r p <C T m T 7 and that the corresponding strength is much smaller than the 
estimates provided by Wiescher et al. First, this level has a known counterpart in 22 Ne at 
£^=6.120 MeV (Endt 1998) with a measured, purely radiative, total width T=29±9 meV, 
that enables to calculate the radiative width of the E x = 5.714 MeV 22 Mg level to be 
T 7 ~23 meV. Since it is found to be greater than the measured total width (16.5±4.4 meV), 
it is likely that r p <r for this level. The second indication comes from experimental data 
on the 21 Ne(d,p) 22 Ne reaction, which provides information on neutron capture on 21 Ne (i.e., 
the mirror counterpart of proton capture on 21 Na). Neutron spectroscopic factors were 
obtained by Neogy et al. (1972) for many 22 Ne levels, from E x =0 to 9.07 MeV, but not 



for the i? x =6.120 MeV one due to its very low population and flat angular distribution. 
Indeed, the direct comparison of peak heights in the experimental spectrum (fig. 1, Neogy 
et al. 1972) shows that, within levels of same spin and transferred angular momentum, the 
cross section for the production of the £^=6.120 MeV level is very small. In addition, the 
corresponding angular distribution (fig. 5, Neogy et al. 1972) is very flat, indicating no 
evidence for a direct component and hence a very small neutron spectroscopic factor. This 
conclusion can be extended to the proton spectroscopic factor in the 22 Na mirror level. The 
6.551 MeV 22 Na level, is assumed (Endt 1990, Endt 1998) to belong to the same isospin 
triplet as the 6.120 MeV 22 Ne level. Its proton spectroscopic factor has been measured by 
Garrett et al. (1971) to be wO.l. This value is not consistent with the Neogy et al. (1972) 
data and would result, for the 5.714 MeV 22 Mg level, in a proton width in excess of the 
measured, 16.5 meV, total width. Hence one may question the isospin triplet assignment 
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made by Endt (1990) and Endt (1998). Moreover, the 5.714 MeV 22 Mg level is withdrawn 
from the triplet in the last paper by Endt (1998). Accordingly, we consider the Wiescher 
et al. estimate for this resonance strength as an upper limit and adopt wy = 2.5, 0.25, 
0.0 meV for upper (r 7 = T p = T/2), recommended value (with the usual 0.1 reduction 
factor) and lower limit. We provide below the updated formula for the rate including the 
uncertainty factors and the direct capture contribution from Caughlan Sz Fowler (1988): 

1.41 x 10 5 T 9 ~ 3/2 exp (-20.739/T 9 1/3 - (T 9 /.366) 2 ) x 

(1. + .02T 9 1/3 + 4.741 T 9 2/3 + 667T 9 + 16.38 T 9 4/3 + 5.858 T 9 5/3 )+ DC (CF88) 

(0 -> 10.) x 40.8. T 9 ~ 3/2 exp (-2.52 /T 9 ) + (E x = 5.714) 

1857. T 9 3/2 exp (-3.95 /T 9 ) + (E x = 5.837) 

(0.1 -»- 10.) x 408. T 9 3/2 exp (-5.49 /T 9 ) (E x = 5.965) 



The third term corresponds to the 5.837 MeV level, assuming that the analog assignment 
made by Endt (1990) is correct. This is however not granted as it is not present anymore 
in Endt (1998). The last term corresponds to the E x =5.965 MeV, J 71 " = + level. Since its 
widths are unknown (save that f 7 <C T Ri T p ), |Wiescher et al. (1986)| assumed a typical 



value for T 7 based on the statistics of Endt (1979). We introduce a factor of ten uncertainty 
to account for the dispersion of the gamma strengths within this statistics, but it has a little 
influence on the rate. The corresponding rates, relative to Caughlan & Fowler (1988) (i.e., 
Wiescher et al. 1986) are shown in Figure 4: the maximum effect occurs between T§ = 0.5 
and 3, well within the nova temperature range. Experimental data of the E x =5.714 MeV, 
j7r_2+ i eve i i s clearly needed to improve the reliability of this rate. 
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A.2. The 22 Mg(p, 7 ) 23 Al and 23 Al(p, 7 ) 24 Si reactions 

As 22 Mg plays a crucial role on the synthesis of 22 Na through its beta decay, its 
destruction by 22 Mg(p,7) 23 Al has to be considered. The contribution of the direct capture 
to the ground state for the 22 Mg(p,7) 23 Al rate has been calculated by [Wiescher et al. 



1986)| . Later on, [Wiescher et al. (1988)| measured the location of the first resonance 
and also calculated its strength (shell model and mirror level). But due to the very low 
Q-value (0.125 MeV), rapid photodisintegration of 23 A1 prevents 22 Mg destruction. Another 
destruction channel has also been proposed: it takes place via two proton captures on 
22 Mg (i.e., proton capture on the small population of 23 A1 through 23 Al(p,7) 24 Si. Gorres, 
Wiescher & Thielemann 1995) and a new rate was proposed by [Herndl et al. (1995)| , based 
on shell model calculations. The location of the first two 24 Si levels has been recently 



determined experimentally (Schatz et al. 1997): it leads to a higher rate than in |Herndl et 



al. (1995"]] due to the lower location of the first resonance (E^ m = 0.141 ±0.031 MeV, instead 



of 320 keV). In our nova models, we have adopted the new limits for the 23 Al(p,7) 24 Si rate 
given by (Bchatz et al. (199T) . 



A.3. The 23 Mg(p, 7 ) 24 Al and 24 Al(p,7) 25 Si reactions 

An estimate to the 23 Mg(p,7) 24 Al rate has been provided by [Wiescher et al. (1986)| 
based on direct capture and three resonances. Since then, Kubono, Kajino & Kato (1995), 
Endt (1998) and Herndl et al. (1998) have discussed analog level assignments, in particular 
with the help of shell model calculations and isobaric multiplet mass equation (Herndl et al. 
1998). For nova nucleosynthesis, only the two first resonances have to be considered. They 
correspond to the E x = 2.349 and 2.534 MeV levels and are assumed to be the analogs of 
the 2.514 MeV, 3+ and 2.563 MeV, 4+(2+) ones in 24 Na (Endt 1998, Herndl et al. 1998). 
In addition to the estimate provided by |Wiescher et al. (1986)| the strength of the second 
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resonance has been calculated (shell model) by |Herndl et al. 199S . However, combining 



statistics of reduced radiative widths (Endt 1979) together with the spectroscopic factor 
reported in [Endt fc Van Per Leun (1978)| , show that its contribution is, in any case, 



negligible with respect to the first one in the domain of temperature considered here. The 
strength of the first resonance is deduced from the analog level (T 7 <C T p ) and hence suffers 
only small uncertainty like the direct capture contribution (see Wiescher et al. 1986). In 
the domain of nova nucleosynthesis, the main uncertainty comes from the determination of 
the energy of the first level: taking E^ m = 0.458 (Kubono et al. 1995), 0.478 (Endt 1998; 
Herndl et al. 1998) or 0.51 MeV (Wiescher et al. 1986) results in a variation of the rate of 
less than a factor of ~ 10 around T 8 = 2.5. To check the importance of this reaction, we 
adopted the highest rate (Kubono et al. 1995) in our test calculations. 

For 24 Al(p,7) 25 Si, we use the rate proposed by [Herndl et al. ( 1995 )| , based on shell 
model calculations. It remains very uncertain, due to the limited experimental spectroscopic 
information on 25 Si. It should have, however, very little effect on Ne-Na leaks, since either 
24 Al(/3+) 24 Mg(p, 7 ) 25 Al or 24 Al(p, 7 ) 25 Si(/3+) 25 Al lead to 25 A1 in the MgAl cycle. 



A.4. The 25 Al(p,7) 26 Si reaction 

25 Al(p,7) 26 Si is important since it leads to the formation of the short-lived isomer 
( 26 Al m ) instead of the long-lived ground state ( 26 A1 5 ). Its rate suffers from large uncertainties 
in the domain of nova nucleosynthesis due mainly to the unknown location of the analog 
of the 26 Mg E x = 6.13 MeV, 3 + level. A level shift similar to the one of its immediate 
neighbor (E x = 6.26 MeV; + in 26 Mg) would bring the 3 + level within the Gamow peak 
with a dramatic influence on the rate. The maximum effect on 26 A1 production is found 
when E^ m (3 + ) = 0.2 ± 0.1 MeV (Coc et al. 1995). New Coulomb displacement energy 
calculations (Iliadis et al. 1996) give E^ m (3 + ) = 0.45±0.1 MeV. However, discrepancies 
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between calculated and experimentally known 26 Si level energies can exceed ±0.1 MeV so 
that the case E^ m (3 + ) = 0.2 cannot be ruled out. (See, e.g., the 23 Al(p,7) 24 Si case above 
where the location of the first resonance is found ~180 keV below the calculated one.) 
Accordingly, for this rate, we used the lower limit, recommended value and higher limit 
provided by Coc et al. (Case A, B and C, respectively). 

A. 5. The 26 Si(p,7) 27 P reaction 

The 26 Si(p,7) 27 P reaction rate has been discussed by |Wiescher et al. (1986)| , and more 
recently by |Herndl et al. (1995)| . Only two excited states are known in 27 P and only one 
(E x = 1.660 MeV, J 71 = 5/2 + ) is low enough to be of astrophysical interest. From the 
comparison with the spectrum of the mirror nucleus, another level, J n = 3/2 + , located at 
E x = 0.985 MeV in 27 Mg, is expected below the 5/2 + one. |Wiescher et al. (1986)| assumed 



a small shift, while shell model calculations (Herndl et al. 1995) led to a much larger shift, 
which resulted in a dramatic increase of the rate by almost four orders of magnitude. The 
calculated strength (Herndl et al. 1995) is in good agreement with the T 7 value (0.8 meV) 
inferred from the lifetime of the E x = 0.985, 27 Mg, conjugate level and with the T p obtained 
from the neutron spectroscopic factor in the conjugate 27 Mg level reported in [Endt fc Van 



Per Leun (1978)| . In consequence, the uncertainty results from the unknown position of 
the resonance. However, it does not affect strongly the rate (a factor of ten at most) for 
1.5 < Tg < 4 down to E^ m ^0.2 MeV due to the compensating effect of the evolution of the 
penetrability Pi p (E r ,s) and exp(—E r /kT) factors with E r small variations. The strength 
of the 5/2 + resonance is known from the |Wiescher et al. (1986)| estimates, from 27 Mg 



spectroscopy and from the |Herndl et al. (1995)| calculations. However, its location (E r 



0.76 MeV) makes its contribution negligible compared to the direct capture one at nova 
temperatures. The latter provides the lower limit for the rate, reflecting the possibility that 
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the J 7 * = 3/2 + resonance energy is lower than E^ m ~0.2 MeV, as considered by |Wicschcr 



et al. (1986)| , while the upper limit is given by the [Herndl et al. (1995)| rate. One should 



note that this effect could be enhanced if the uncertain 25 Al(p,7) 26 Si rate is close to its 
maximum value (see above). 

A.6. The 27 P(p, 7 ) 28 S, and 27 Si(p, 7 ) 28 P reactions 

Only the ground state of 28 S is known but from the level scheme of the conjugate 
nucleus, 28 Mg, one can infer that no resonance is expected for E^ m ^1.4 MeV. Hence, at 
the temperatures considered in the text, the 27 P(p,7) 28 S reaction proceeds only through 
direct capture. The corresponding rate has been obtained by [Herndl et al. ( 1 995 )| from shell 
model calculations. It is about a factor of ten higher than the Caughlan & Fowler (1988) 
one and can be considered as reliable enough (i.e. within an order of magnitude), due to 
the expected absence of resonances. 

The 27 Si(p,7) 28 P reaction rate given in Caughlan & Fowler (1998) comes from the 
analysis of IWicscher et al. (1986) . They considered five resonances corresponding to 



the 28 P levels between E x = 2.143 and 2.628 MeV, out of which, for three levels, the 
spins are uncertain. From the assumed conjugate levels in 28 A1, they extracted radiative 
widths and proton widths using spectroscopic factors obtained by neutron transfer on 
27 A1. More recently, the good agreement between calculated (Endt & Booten 1993) and 
experimental (Endt 1990) level schemes, spectroscopic factors and radiative widths, for 
A = 28 nuclei, gives confidence to the 28 A1 and 28 P conjugate level assignment (see also 
lliadis et al. 1998| ). The first level above proton threshold, E x = 2.104 MeV, J w = 2 + , 



was not considered by [Wiescher et al. (1986)| , but its estimated proton width is however 



too small (r p ^ rjf ~ 10 20 eV) to contribute to the rate for T 8 ^ 0.3. The uncertainty 
remains limited to transposition of 28 A1 data to 28 P. 
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Table 1. Initial composition of the envelope (up to Si), assuming 50% mixing with the 

ONe white dwarf core 



Nuclei Mass fraction 



iff 


3.5E-1 


3 He 


1.5E-5 


4 He 


1.4E-1 


6 Li 


3.2E-10 


7 Li 


4.7E-9 


9 Be 


8.3E-11 


io B 


5.3E-10 


n B 


2.4E-9 


12 C 


6.1E-3 


13 C 


1.8E-5 


14 N 


5.5E-4 


15 N 


2.2E-6 


16Q 


2.6E-1 


17 


1.9E-6 


18Q 


1.1E-5 


19p 


2.0E-7 


20 Ne 


1.6E-1 


21 Ne 


3.0E-3 


22 Ne 


2.2E-3 


23 Na 


3.2E-2 


24 Mg 


2.8E-2 


25 Mg 


7.9E-3 


26 Mg 


5.0E-3 
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Table 1 — Continued 

Nuclei Mass fraction 

27 A1 5.4E-3 

28 Si 3.3E-4 
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Table 2. Test models 



Nuclear reaction 


Old rate 


Test rate 


X( 22 Na) tost 
X( 22 Na) old 


X( 26 Al) tcst 

X( 2 6Al) ld 


1.15 M ONe 


21 Na(p, 7) 22 Mg 


CF88 


CF88/100 


2.3 


1.1 


25 Al(p, 7 ) 26 Si 


Wie86 


Coc95, case A 


1 


1.2 




Wie86 


Coc95, case C 


1 


0.5 


23 Na(p, 7 ) 24 Mg 


CF88 


CF88+GWR89 


1.2 


1.3 


26 AP(p, 7 ) 27 Si 


Vog89 


Vog89 
l/3xres(0.188 MeV) 


1 


1.9 


26 Al m (p, 7 ) 27 Si 


CF88 


CF88xl00 


1 


1 


26 Si(p, 7 ) 27 P 


Wie86 


Her95 


1 


1 


26 Mg(p, 7 ) 27 Al 


11190 


Ili90+Cha90 


1 


1 


1.25 M ONe 


21 Na(p, 7 ) 22 Mg 


CF88 


CF88/100 


3 


1.1 


22 Na(p, 7 ) 23 Mg 


CF88 


Ste96 


3 


1.2 


23 Mg(p, 7 ) 24 Al 


Wie86 


KKK95 


1.1 


1.1 


1.35 M ONe 


21 Na(p, 7 ) 22 Mg 


CF88 


CF88/100 


1.2 


1 


23 Al( P)7 ) 24 Si 


Wor94 


Sch97 


1 


1 



References to the reaction rates: Wie86 (Wiescher et al. 1986), CF88 
(Caughlan & Fowler 1988), GWR89 (Gorres, Wiescher & Thielemann 1989), 
Vog89 (Vogelaar 1989), Cha90 (Champagne et al. 1990), Ili90 (Iliadis et 
al. 1990), Wor94 (van Wormer et al. 1994), Coc95 (Coc et al. 1995), 
Her95 (Herndl et al. 1995), KKK95 (Kubono, Kajino & Kato 1995), Ste96 
(Stegmiiller et al. 1996), Sch97 (Schatz et al. 1997) 
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Table 3. Reaction rates adopted for the calculation of Maximum (A), Recommended (B), 
and Minimum (C) 22 Na- 26 Al production compared with those adopted in Jose & Hernanz 

(1998) 



Reaction 



JH98 



A 



B 



C 



21 T 



Na(p, 7 ) 22 Mg CF88 



This work 
Lower limit 



This work 
Recommended 



This work 
Upper limit 



21 Ne(p, 7 ) 22 Na 


CF88 


CF88 


CF88 


CF88 




with f=0.1 


f=0.01 


f=0.001 


f=0 


22 Na(p, 7 ) 23 Mg 


CF88 


Ste96 


NACRE 


Ste96 






Lower limit 


Recommended 


Upper limit 


25 Al(p, 7 ) 26 Si 


Wie86 


Coc95 


Coc95 


Coc95 






Case A (Lower) 


Case B (Recommended) 


Case C (Upper) 


23 Na(p, 7 ) 24 Mg 


CF88 


CF88+GWR89 


CF88+GWR89 


CF88+GWR89 






with f=l 


f=0.1 


f=0 


23 Mg(p, 7 ) 24 Al 


Wie86 


KKK95 


KKK95 


KKK95 


26 Mg(p, 7 ) 27 Al 


Ili90 


Cha90+Ili90 


Cha90+Ili90 


Cha90+Ili90 






with f=l 


f=0.1 


f=0 


26 Al m (p, 7 ) 27 Si 


CF88 


CF88 


CF88 


CF88 


27 Al(p,a) 24 Mg 


Cha88 


Cha88+Tim88 


Cha88+Tim88 


Cha88+Tim88 






with f=l 


f=0.1 


f=0 


23 Al(p, 7 ) 24 Si 


WW80 


Sch97 


Sch97 


Sch97 






Lower limit 


Recommended 


Upper limit 


26 Si( P , 7 ) 27 P 


Wic86 


Her95 


Hcr95 


Her95 


26 A19(p, 7 ) 27 Si 


Vog89 


Coc95 


Coc95 


Coc95 






with fl=f2=0 


fl=f2=0.1 


fl=f2=l 



Additional references to the reaction rates: WW80 (Wallace & Woosley 1980), Cha88 
(Champagne et al. 1988), Tim88 (Timmcrmann ct al. 1988), NACRE (Angulo ct al. 1998). f, 
fl and f2 represent uncertainty factors in the analytic rates. 
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Table 4. Main characteristics of the 1.15 M ONe nova models compared with model 

ONe3 (Jose & Hernanz 1998) 

ONe3 ONell5A ONell5B ONell5C 



Reaction Network 


Old 


A 


B 


C 


AM env (1(T 5 


M ) 


3.2 


3.2 


3.2 


3.2 


T CIO 8 

-'-max V w 


K) 


2.19 


2.33 


2.31 


2.30 


AM ejec (1(T 5 


■M ) 


1.9 


2.6 


2.6 


2.6 


K (10 45 erg) 


1.2 


1.5 


1.5 


1.5 


20 Ne 




1.8E-1 


1.6E-1 


1.7E-1 


1.8E-1 


21 Ne 




3.0E-5 


2.0E-4 


8.7E-5 


2.5E-5 


22 Ne 




1.7E-3 


1.6E-3 


1.5E-3 


1.5E-3 


22 Na 




5.3E-5 


5.3E-4 


2.7E-4 


1.2E-4 


23 Na 




7.5E-4 


1.3E-3 


1.4E-3 


1.4E-3 


24 Mg 




1.0E-4 


9.6E-4 


4.1E-4 


3.5E-4 


25 Mg 




2.9E-3 


6.1E-3 


3.9E-3 


1.8E-3 


26 Mg 




3.4E-4 


8.3E-4 


6.4E-4 


1.2E-3 


26 Al 




9.3E-4 


1.2E-3 


8.0E-4 


3.3E-4 


27 A1 




4.5E-3 


4.8E-3 


3.5E-3 


3.4E-3 


28 Si 




5.4E-2 


6.8E-2 


5.8E-2 


5.6E-2 
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Table 5. Main characteristics of the 1.25 M ONe nova models compared with model 

ONe5 (Jose & Hernanz 1998) 

ONe5 ONel25A ONel25B ONel25C 



Reaction Network 


Old 


A 


B 


C 


AM env (1(T 5 


M ) 


2.2 


2.2 


2.2 


2.2 


T CIO 8 

-'-max V w 


K) 


2.44 


2.54 


2.51 


2.51 


AM ejec (1(T 5 


■M ) 


1.4 


1.8 


1.8 


1.8 


K (10 45 erg) 


1.4 


1.5 


1.5 


1.5 


20 Ne 




1.8E-1 


1.6E-1 


1.7E-1 


1.7E-1 


21 Ne 




3.5E-5 


3.7E-4 


1.4E-4 


3.0E-5 


22 Ne 




1.0E-3 


9.8E-4 


9.1E-4 


8.6E-4 


22 Na 




9.6E-5 


6.9E-4 


3.5E-4 


1.9E-4 


23 Na 




1.4E-3 


1.6E-3 


1.9E-3 


1.9E-3 


24 Mg 




2.0E-4 


9.0E-4 


4.1E-4 


3.9E-4 


25 Mg 




2.4E-3 


4.6E-3 


3.2E-3 


1.3E-3 


26 Mg 




2.8E-4 


5.9E-4 


5.1E-4 


1.1E-3 


26 Al 




5.4E-4 


8.8E-4 


6.5E-4 


2.0E-4 


27 A1 




2.0E-3 


2.9E-3 


2.2E-3 


2.2E-3 


28 Si 




5.6E-2 


6.8E-2 


5.8E-2 


5.6E-2 
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Table 6. Main characteristics of the 1.35 M ONe nova models compared with model 

ONe6 (Jose & Hernanz 1998) 

ONe6 ONel35A ONel35B ONel35C 



Reaction Network 


Old 


A 


B 


C 


AM env (10- 5 


M ) 


0.54 


0.54 


0.54 


0.54 


T CIO 8 

-'-max V w 


K) 


3.24 


3.22 


3.31 


3.23 


AM ejec (10- 5 


■M ) 


0.44 


0.44 


0.44 


0.44 


K (10 45 erg) 


0.9 


1.0 


1.0 


1.0 


20 Ne 




1.5E-1 


1.3E-1 
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Fig. 1. — Snapshots of the evolution of 20 ' 21 ' 22 Ne, 21 ' 22 ' 23 Na, and 22,23 Mg (in mass fractions) 
along the accreted envelope for a 1.25 M ONe novae accreting at a rate M = 2 xlO~ 10 
M Q .yr _1 . The mass coordinate represents the mass below the surface, relative to the total 
mass. The successive panels, from top to bottom, correspond to the time for which the 
temperature at the burning shell reaches 5 xlO 7 , 7 xlO 7 , 10 8 , 2 xlO 8 , T max = 2.44 
xlO 8 , plus 3 panels corresponding to the last phases of the explosion, when the white dwarf 
envelope has already expanded to a size of R w d ~ 10 9 , 10 10 and 10 12 cm, respectively. The 
arrow indicates the base of the ejected shells. 

Fig. 2. — Main nuclear paths in both NeNa and MgAl cycles. 

Fig. 3.— Same as Fig. 1, for the evolution of 24 > 2 5,26 Mg) 26 M g,m^ 25,27^ and 26,27 Si 

Fig. 4. — Comparison between the nuclear reaction rate for 21 Na(p,7) 22 Mg provided in this 
work, with the analytic rate given in Caughlan & Fowler (1988). 
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